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Abstract 

The LULU operators, well known in the nonlinear multiresolution 
analysis of sequences, are extended to functions defined on continuous 
domain, namely, a real interval ft C K. Similar to their discrete coun¬ 
terparts, for a given <5 > 0 the operators Lg and Us form a fully ordered 
semi-group of four elements. It is shown that the compositions Ls o Us 
and Us o Ls provide locally d-monotone approximations for the bounded 
real functions defined on ft. The error of approximation is estimated in 
terms of the modulus of nonmonotonicity. 


1 Introduction 

The LULU operators remove impulsive noise before a signal extraction from 
a sequence. They are computationally convenient and conceptually simpler 
compared to the the median smoothers usually considered to be the ’’basic” 
smoothers. The LULU operators have particular properties, e.g. they are fully 
trend preserving, | 3 ], preserve the total variation, [ 2 ], etc., which make them 
an essential tool for multiresolution analysis of sequences. Furthermore, it was 
demonstrated during the last decade or so that these operators, being specific 
cases of morphological filters, [^, have a critical role in the analysis and com¬ 
parison of nonlinear smoothers, @]. 

We extend the LULU theory from sequences to functions on a continuous 
domain, namely, a real interval Ul. The existing LULU theory can be considered 
as a particular case in this new development, since the discrete LULU operators 
can be equivalently formulated for splines of order 1 or order 2 on the integer 
partition of real line, the one-to-one mapping being given by the B-spline basis. 

Given a sequence ^ and n G N the operators and Un are defined 

as follows 

{LnOi = max{min{^j_„, ...,^i}, ...,min{^i, i € N 

{UnOi = min{max{^j_„,..., ..., maxj^i,..., * € N 
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In analogy with the above discrete LULU operators, for a given ^ > 0 the ba¬ 
sic smoothers Ls and Us in the LULU theory are defined for functions on fl 
through the concepts of the so called lower and upper d-envelopes of these func¬ 
tions. These definitions are given in Section 2, where it is also shown that the 
operators Ls and Us preserve essential properties of their discrete counterparts. 
In particular, the operators Ls and Us generate through composition a fully 
ordered four element semi-group also called a strong LULU structure. This 
issue is dealt with in Section 3. In section 4 we define the concept of local 5- 
monotonicity and show that the compositions Ls o Us and Us o Ls are smoothers 
in the sense that the resulting functions are locally (5-monotone. The errors of 
approximation of real functions / by the these compositions are estimated in 
terms of the modulus of nonmonotonicity ^{f,S). 


2 The basic smoothers Ls and Us 

Let ^(U) denote the set of all bounded real functions defined on the real interval 
U C R. Let Bs{x) denote the closed ^-neighborhood of x in U, that is, Bs{x) = 
{y : \x — y\ < (5}. The pair of mappings J, S : ^(U) ^ ^(U) defined by 

I{f){x) = snpini{f{y) :y G Bsix)}, X en, (1) 

(5>0 

= inf sup{f{y) : y G Bs{x)}, x gLI. (2) 

(5>0 

are called lower Baire, and upper Baire operators, respectively, |S]. We consider 
on ^(U) the point-wise defined partial order, that is, for any f,gG yI(U) 

f < g f{x) < g{x), X Gn. (3) 

Then the lower and upper Baire operators can be defined in the following 
equivalent way. For every / G ■4(U) the function /(/) is the maximal lower 
semi-continuous function which is not greater than /. Hence, it is also called 
lower semi-continuous envelope. In a similar way, S{f) is the smallest upper 
semi-continuous function which is not less than / and is called the upper semi- 
continuous envelope of /. In analogy with /(/) and S{f) we call the functions 

Is{f){x) = inf{/(y) : y G Bs{x)}, x G n, (4) 

Ssif)ix) = sup{/(y) : y G Bs{x)}, x G Ll. (5) 

a lower ^-envelope of / and an upper (5-envelope of /, respectively. 

It is easy to see from Q and © that for every 61,62 >0 

L\ ^ ^82 ^S\-\-S 2 5 O Us 2 *^i5i-t-{52 (^) 

Furthermore, the operators Is and 5^, ^ > 0, as well as I and S are all monotone 
increasing with respect to the order ®, that is, for every f,gG .4(U) 

f<g=^ h{f) < Is{g), Ssif) < Ssig), I{f) < I{,g), S(,f) < S{g). (7) 
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The following operators can be considered as continuous analogues of the 
discrete LULU operators given in the Introduction: 

Lg = Ss o Is , Ug = Is_ o S± . 

2 2 2 2 

We will show that these operators have similar properties to their discrete coun¬ 
terparts. Let us note that they inherit monotonicity with respect of the func¬ 
tional argument from the operators Ig and Sg, o, that is, for f,gG 

f<g^ Lsif) < Ls{g), Ug{f) < Ug{g). ( 8 ) 

Theorem 1 For every fGA{Ll) and S>0 we have Lg{f) < /, Ug{f) > f. 

Proof. Let / G .4(n), 5 > 0. For any a; G it follows from the definition of 
Ig that U{f){y) < f{x), y G Bs{x). Therefore Lg{f){x) = Ss{Is{f)){x) = 
sup{Is{f){y) : y G Bs{x)} < f{x), x £ Ll. The second inequality in the theorem 
is proved in a similar way. ■ 

Theorem 2 The operator Lg is monotone increasing on 5 while the operator 
Ug is monotone decreasing on 5, that is, for any f G and 0 < Si < 62 we 

have Lg^{f) < Lg^{f), Ug,{f) > Ug^{f). 

Proof. Let (52 > > 0. Using the properties the operator Lg^ can be 

represented in the form Lg^ = Ss^ o Is^ = Ss-i o S'sa-^i o Fz-h ° = Ss^ o 

2 2 2 2 2 2 2 
Lg^-g^ o I^. It follows from Theorem ^ that for every / G A{Ll) we have 
2 

Lg2-gi{Inif)) < Hence using the monotonicity of the operator Sg 

given in Q we obtain Lg^{f) = 5^(/))) < Snilnif)) = Lg^f), 

2 2 2 2 

/ G A{Ll). The inequality Ug^{f) > Ug^{f) is proved in a similar way. ■ 

The next lemma is useful in dealing with compositions of Ig and Sg. 

Lemma 1 We have Ig o Sg o Ig = Ig, Sg o Ig o Sg = Sg. 

Proof. Using the monotonicity of Ig, see o, and Theorem □ for / G ^(U) 
we have {Ig o Sg o Ig){f) = Ig{L2g{f)) < h{f)- On the other side, applying 
Theoremnto U2g we obtain {Ig o Sg o Ig){f) = U2g{Ig{f)) > Is{f)- Therefore 
{Ig o Sg o Ig){f) = Ig{f), f G .4(f2). The second equality is proved similarly. ■ 

Theorem 3 For every 61,62 > 0 we have Lg^oLg.^ = Lmax{i5i,i52} ^xid = 

Umax{(5i,(52} ■ 

Proof. We will only prove the hrst equality since the proof of the second one 
is done in a similar manner. Let first ^2 > ^1 > 0. Using the property and 
Lemma [n we obtain 

o {S oJsi^o(^Ss2 ^ IS2 ^ {S o J o S Si ^ o ^ S S2 — o T S2 ^ 

~T ~T ~T ~T ~T ~T ~T 2 “ 

— S 5;]^ O S <52 —5]^ O 1 82 — S 82 O <52 — ■ 

“ 2 “ ~ ~T 
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If i5i > (52 > 0 in a similar way we have 

® -^(52 (*^(51 ® ^Si') ® {^62 ® ^ 52 ) (*^5i ^61 — 82 ^ ^ (-^^2 ^ *^52 ® ^ 82 ^ 

= Ssi o Is^_ 82 o I82 = Ssi o Is^ = Ls^. 

The proof in the case when ^2 = ^r > 0 follows from either of the above identities 
where SS 2-81 or Isi -82 respectively are replaced by the identity operator. ■ 
Important properties of smoothing operators are their idempotence and co- 
idempotence. Hence the significance of the next theorem. 

Theorem 4 The operators Lg and Ug are both idempotent and co-idempotent, 
that is, Lgo Lg = Lg, Ug oUg = Ug, {id — Lg) o {id — Lg) = id — Lg, {id — Ug) o 
{id — Ug) = id — Ug, where id denotes the identity operator. 

Proof. The idempotence of Lg and Ug follows directly from Theorem |31 The 
co-idempotence of the operator Lg is equivalent to Lg o {id — Lg) = 0. Using the 
first inequality in Theorem ^ one can easily obtain Lg o {id — Lg) > 0. Hence, 
for the co-idempotence of Lg it remains to show that Lg o {id — Lg) < 0. Assume 
the opposite. Namely, there exists a function / G A(H) and x G Lt such that 
{Lg o {id — Lg)){f){x) > 0. Let e > 0 be such that {Lg o {id — Lg)){f){x) > 
e > 0. Using the definition of Lg the above inequality implies that there exists 
y G Bs{x) such that for every z G Bs{y) we have {id — Lg){f){z) > e, or 
equivalently 

f{z) > Lg{f){z) -be, zG Bs {y). (9) 

For every z G B^{y) we also have Lg{f){z) > Is{f){y) = inf{/(t) : t G B^{y)}. 
Hence there exists t G B^{y) such that f{t) < Is{f){y) -b e < Lg{f){z) + e, 
-2 G Bs{y). Taking z = t in the above inequality we obtain f{t) < Lg{f){t) -be, 
which contradicts ©■ The co-idempotence of Ug is proved in a similar way. ■ 

Example 1 The figures below illustrate graphically the smoothing effect of the 
operators Lg, Ug and their compositions. The graph of function f is given by 
dotted lines. 




The functions Lg{f) and Ug{f) 
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The functions {Ls o Us){f) and {Us o Ls){f) 

The operator Ls smoothes the function / from above by removing sharp 
picks while the operator Us smoothes the function / from below by removing 
deep depressions. The smoothing effect of the compositions Ls o Us and Us o Ls 
can be described in terms of the local i5-monotonicity discussed in the Section 4. 

Note that LsoUs and UsoLs resolve ambiguities in a different way; LsoUs treats 
oscillations of length less then 6 as upward impulses and removes them while 
Us o Ls considers such oscillations as downward impulses which are accordingly 
removed. The inequality {Us o Ls){f) < {Ls o Us){f) which is observed here will 
be proved in the next section for / £ y4(n). 

3 The LULU semi-group 

In this section we consider the set of the operators Ls and Us and their compo¬ 
sitions. For operators on y4(n) we consider the point-wise defined partial order. 
Namely, for operators P, Q on ^(17) we have 

P<Q ^ P{f) < Q(/), / £ ^(f7). 

Then the inequalities in Theorem ^ can be represented in the form 

Ls ^ id < Us, (10) 

where id denotes the identity operator on ^(17). 

Theorem 5 For any S > 0 we have Us o Ls < Ls o Us- 

Proof. Let / £ and let x £ 17. Denotep = {LsoUs){f){x) = Ss {Is{S^{f))){x). 
Let e be an arbitrary positive. For every y £ Bs^{x) we have 

Is{S^{f)){y) <p <p +e. (11) 

Case 1. There exists z £ Bs{x) such that Ss{f){z) < p + e. Then f{t) < 
p + e for t £ Si (z), which implies that/£(/)(t) < p-|-£ for t £ ^^(z). Hence 
Ss{P^{mz) < p + e. Then {Us o Ls){f){t) = n{Ss{Is{m{t) < p + e for 
t £ Bs (z). Since x £ Sa (z), see the case assumption, from the above inequality 
we have {Us o Ls){f){x) < p + e. 
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Case 2. For every z G Bs (x) we have S6{f){z) > p+e. Denote D = |z G D : Ss{f){z) < p + e 
We will show that for every z G Bs(x) we have 

Bs{z)r\D^% (12) 

Due to the inequality m we have that 03 holds for every z G Bs{x). Let 
z G Bs{x) and let z > a; + |. This implies that a; + | G Using the inequality 
itTTll for 2 / = a; + I as well as the case assumption we obtain that the set 
(a; + |,a; + ^] n D is not empty. Then Bs{z) H D D (a; + |,a:+^] H D ^ 0. 

For z < X — I condition 03 is proved in a similar way. Hence 03 holds for all 
z G Bs{x). Let z G Bs{x) and v G Bs{y)r\D. Since v G D we have f{t) < p + e, 
for t G Bs (v). Using that Bs{z) n Bs{v) 0 we obtain that Is{f){z) <p + e, 
z G Bs{x). Therefore Ss{Is{f j){x) <p + e. Then 

(Us o Ls){f)ix) = Is (Ssils (mix) < Ss{Is (/))(x) <p + e. 

Combining the results of Case 1 and Case 2 we have {Us o Ls){f){x) < p + e. 

Since e is arbitrary this implies that {Us o Ls){f){x) <p= {Ls o Us){f){x). ■ 

Theorem 6 For a given <5 > 0 the operators Ls o Us and Us o Ls are both 
idempotent. 

The proof is an immediate application of Lemma 1. 

Theorem 7 We have Us o Ls o Us = Ls o Us, Ls o Us o Ls = Us o Ls, (5 > 0. 

Proof. Using the inequalities 03 and the monotonicity of the operators Ls, 

Us, see we obtain Us o Ls o Us > id o Ls o Us = Ls o Us- For the proof of 
the inverse inequality we use Theorem 13 and the idempotence of Us as follows: 

UsoLsoUs = {Us o Ls) oUs< {Ls o Us) oUs = Lso {Us o Us) = LsoUs 

Therefore Us o Ls o Us = Ls o Us- The second equality is proved in a similar 
way. ■ 

It follows from Theorems 0 and [3 that for a fixed 5 > 0 every composition 
involving finite number of the operators Ls and Us is an element of the set 
{Ls, Us, Us o Ls,Ls o Us}- Hence the operators Ls and Us form a semi-group 
with a composition table as follows: 



Ls 

Us 

Us 0 Ls 

Ls 0 Us 

Ls 

Ls 

Ls 0 Us 

Us 0 Ls 

Ls 0 Us 

Us 

Us 0 Ls 

Us 

Us 0 Ls 

Ls 0 Us 

Us 0 Ls 

Us 0 Ls 

Ls 0 Us 

Us 0 Ls 

Ls 0 Us 

Ls 0 Us 

Us 0 Ls 

Ls 0 Us 

Us 0 Ls 

Ls 0 Us 


Furthermore, an easy application of Theorem 0 shows that this semi-group is 
completely ordered. Namely, we have Ls < Us o Ls < Ls o Us < Us- 
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4 Locally (^-monotone approximations 

Definition 1 Let 5 > 0 and a function f G A(f2) be given. 

(i) The funetion f is called upwards 5-monotone if for every interval [x,y] C 
Vi with y — X < 5 we have sup f{z) = m.ax{f{x), f{y)}. 

ze[x,v] 

(a) The function f is called downwards 5-monotone if for every interval 
[x,y\ C LI with y — x <5 we have inf f{z) = min{/(a;),/(y)}. 

ze[x,y] 

(Hi) The function f is called locally 5-monotone if it is both downwards 5- 
monotone and upwards 5-monotone. 

The name locally J-monotone reflects the following characterization: 

/ is locally _ On any interval [x^y] CO, y — x<5, f is either , , 

^-monotone monotone increasing or monotone decreasing ^ 

Theorem 8 For every i5 > 0 and f C ^(O) the funetion Ss^f) is upwards 
5-monotone while the funetion Is{f) is downwards 5-monotone. 

Proof. Let 5 > 0, f G .4(0) and [x,y] Q Ll, y — x < 5. Denote g = Ss{f). It is 
easy to see that for every z C [x,y] we have Bs{z) C B^{x) Ui?|(y). Therefore 

g{z) = snp{f (t) :t G B6 {z)} < sup{f (t) : t G B^{x) U Bs{y)} 

= max{sup{/(t) : t G (x)},sup{/(t) : t G B^{z)}} 

= max{g{x),g{y)}, 

which shows that function g is upwards 5-monotone. The downwards 5-monotonicity 
of Us{f) is proved in a similar way. ■ 

Corollary 1 For every 5 > 0 and f G A(fl) the function Ls{f) is upwards 
5-monotone while the funetion Us{f) is downwards 5-monotone. 

Theorem 9 For every 5 > 0 and f C A{Ll) the functions {Us o Ls){f) and 
{Ls o Us){f) are both locally 5-monotone. 

The proof follows from Corollary ^ and the composition table in the preceding 
section. Theorem El shows that the operators Us o Ls and Ls o Us provide 
locally 5-monotone approximations to the functions in ^(f2). The error of the 
approximation can be estimated in terms of the modulus of nonmonotonicity. 

Let us recall the definition. 

Definition 2 Let fGA{Ll). The mapping fj.{f, •) : R’*'—i-R+UjO} given by 

= ^ sup sup (|/(xi)-/(x)H/(a: 2 )-/(a:)|-|/(xi)-/(x 2 )|) 

Xi ^2 & LI X G [xi, X 2 ] 

0<X2—X2 <5 

is called modulus of nonmonotonicity of f. 
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The locally J-monotone functions can be conveniently characterized through the 
modulus of nonmonotonicity. For any / G ■4(0) and 5 > 0 we have 

/ is locally (5-monotone /r(/, (5) = 0 (14) 

We will derive error estimates first in the case when fl = R. It will prove 
useful to consider the upper semi-continuous envelope of the modulus of non¬ 
monotonicity. Let / G w4(n). Using that is monotone increasing with 

respect to S the upper semi-continuous envelope of /i can be represented as 
A(/,<5) = = lim£^o+ + 

Theorem 10 Let f G yl(M) and (5 > 0. Then f{x) — Ls{f){x) < jj,{f,6), 
Us{f)ix) - f{x) < /!(/, S), X gLI. 

Proof. Let a: £ M. Denote p = Ls{f){x). li p = f{x) the first inequality of the 
theorem holds. Assume that p < f{x). Let 77 > 0 be such that p + rj < f{x). 
Then we have 

I^if){y) <P + V, y G B6 {x). (15) 

Denote Di = {z < x : f{z) < p + rj}, D 2 = {z > x : f{z) < p + p}, zi = supDi, 

Z 2 = swpD 2 . Using the inequality 11511 with y = x — ^ and y = x + ^ we obtain 
that Di and respectively D 2 are not empty and that x — S < zi < x < Z 2 < x+S. 
Therefore Z 3 = G 2 Lh|± 5 ] ^ Bsix). Then the inequality ^ 

implies that Bs n {Di U D2) 0 . Hence Z2 — Z3 = Z3 — zi < |. 

Let e > 0 be arbitrary. The neighborhood B 6 +e(z 3 ) has nonempty intersec¬ 
tions with both Di and ZI 2 . Let ti G B^ n Di and t 2 G B^+e fl D 2 . We have 

2 2 

t 2 —ti < S + 6 and X G [fi, 12 ] ■ From the definition of the modulus of nonmono¬ 
tonicity we have If(ti)-f(x)l + lf(t 2 )-f(x)i-lf(ti)-f(t 2 )l < 2p(f,S + e). On 
the other side |/(ti) - /(x)| -f 1/(0) - /(a;)| - |/(ti) - /( 0 )| = 2 f(x) - /(ti) - 
/(0)-|/(0)-/(0)| = 2/(x)-max{/(ti),/(0)} > 2f(x)-2(p + p). Therefore 
f(x) — p — p < y,{f, (5 -I- e). Going with e to 0 we obtain f{x) —p< 6) -I- p. 

Since p is arbitrary small this implies the first inequality of the Theorem. The 
second inequality is proved in a similar manner. ■ 

Using Theorem ^1 as well as la and Corollary ^ we have the following 
characterization of the hxed points of operators Ls and Us- For any / G A(M) 

Kf,S)=0 iLs{f)=f, Usif) = f) y{f,S)=0 

Theorem 11 LetS > 0 and f G A(R) . Then fi{Ls{f), 6) < fi{f,d), p{Ls{f),6) < 
Kf,s), p{Usif),d)<pif,s), piUsif),s)<Kf,s). 

Proof. We will prove the inequalities for Ls{f) since the ones for Us{f) are 
proved in a similar way. Denote g = Ls{f). Let [xi, 0 : 2 ] be an arbitrary interval 
of length at most 5 and let x G [xi,X 2 ]- We consider the number q = \ f{xi) — 
/(x)| -I- |/(x 2 ) — f{x)\ — |/(xi) — /(x 2 )|. According to Corollarynthe function 
g is upper 6 monotone, which implies that g(x) < uiax{g{xi), g{x 2 )}- If we 
also have g{x) > min{g(xi), ( 7 (x 2 )}, then the number q is zero and the first 
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inequality of the Theorem is trivially satisfies. Let g(x) < min{g(a:i), (?(x 2 )} 
and let ry > 0 be such that g{x) + ij < min{(/(xi), ( 7 ( 3 : 2 )}. The number q can 
then be represented in the form q = 2 (min{( 7 (xi), 5 ( 0 : 2 )} — g{x)). If we assume 
that f{y) > g{x)+r] for all y £ Bs (x), then g{x) > Is {f){x) > g{x)+r], which is 
a contradiction. Therefore, there exists y £ Bs{x) such that f{y) < g{x) + 5 . If 
y < xi then using that Bs (a:i) C Bs {y) U Bs (x) we obtain B (z) < g{x) + g for 
all 2 : G Bs{xi) which implies 5 ( 0 : 1 ) = Ss{Is{f)){xi) < g{x) + g < 5 ( 0 : 1 ). This 
contradiction shows that y > Xi. In a similar way we show that y < X 2 - Using 
also the first inequality of Theoremn]we have q = 2 (min{ 5 (a:i), 5 ( 3 : 2 )} — 5 ( 0 :)) < 
2(min{/(o:i), 7 ( 0 : 2 )} — /(y) — 2g) < fi{f, 6 ) — 2g. Using that 5 can be arbitrary 
small we obtain q < fi{ f, S). Since the interval [ 0 : 1 , 0 : 2 ] of length at most S and 
X £ [xi,o: 2 ] are arbitrary this implies that g{g,S) < yL{f, 6 ). The inequality 
^) < A(/i 5) follows from the monotonicity of the operator S, see 0 . ■ 

In the next theorem we give estimates for the error of approximation of a 
function / £ ^(R) in terms of the supremum norm denoted here by 11 • 11. 

Theorem 12 Let J > 0 and f £ ^(R) . Then \\f — {Ls o Us){f)\\ < 
\\f-{UsoLs){f)\\<fi{f,S). 

Proof. Applying Theorems Uni and ITTI we obtain (Ls oUs){f){x) < Us{f){x) < 
fix) + fiif, 6 ) a,iid {Ls oUs){f)ix) >Usif){x) - fl{Lsif),S) > f {x) - fiif, 6 ), 
which implies the first inequality of the Theorem. The second inequality is 
proved in a similar way. ■ 

Error estimates similar to Theorem II 2l ca,n be derived in case of U being finite 
or semi-finite interval using a modification of the modulus of nonmonotonicity. 
For simplicity we will only consider the case Tl = [a,h], a,h £ R. 

Definition 3 Let /gA(U). The mapping A(/i') ^ R'''^R'''U{0} given by 

fL{f,6) = sup < gif, 6), sup (|/(xi)-/(x 2 )|, sup (|/(xi)-/(x 2 ) 

} a:i,2e[o,o-|-|] a:i,2e[fc—|.f)] 

is called modified modulus of nonmonotonicity of f. 

This modulus is similar to the corrected modulus of nonmotonicity in ^ . 

Theorem 13 Let f £ A(R) and <5 > 0 . Then fix) — Lsif)ix) < gif, 6 ), 
Usif)ix) - fix) < gif, S), X £fl. 

The proof is similar to the proof of Theorem llOl 

It is easy to see that for any / G Ala, b] the functions Lsif) and Usif) are 
constants on each of the intervals [a, a -I- |] and [b — |]. Therefore, using also 
Theorem im we have giLsif), 6 ) = giLsif),S) < gif, 6 ) < gif,5). In the same 
way we obtain giUsif),S) < gif,5). Hence the modified modulus satisfies the 
similar inequalities to the ones given in Theorem 1111 for g and g. Using these 
inequalities and Theorem ll3l we obtain the error estimates in the next theorem, 
which are similar to the ones in Theorem [El 

Theorem 14 Let S > 0 and f £ A(R) . Then \\f — (L^ o Us)if)\\ < TifA)j 
\\f-iUsoLs)if)\\<gif,5). 
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5 Conclusion 


In this paper we extended the LULU operators from sequences to real functions 
defined on a real interval using the lower and upper (5-envelopes of functions. 
The obtained structure, although more general than the well known LULU 
structure of the discrete operators, retains some of its essential properties. For 
a fixed <5 > 0 the compositions Ls o Us and Us o Ls provide locally (5-monotone 
approximations for real functions, the error of approximation being estimated 
in terms of the modulus of nonmonotonicity of the functions. Further properties 
of the LULU operators for functions on continuous domains, e.g. trend preser¬ 
vation, will be investigated in the future. Generalizing the theory to functions 
on multidimensional domains is still an open problem. 
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